Ground state properties of a fully frustrated quantum spin system 
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We find that ground states of the quantum Heisenberg antiferromagnet on the geometrically frus- 
trated pyrochlore checkerboard lattice are singlets and can be expressed in terms of positive matrices. 
The magnetization at zero external field vanishes for each frustrated tetrahedral unit separately and 
there is an upper bound of 1/8 in natural units on the susceptibility both for the ground state and 
at finite temperature. These results are the first exact ones in this field and generalize to some other 
lattices; the approach is also of interest for other spin systems. 
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Geometrically frustrated spin systems are known to 
have many interesting properties that are quite unlike 
those of conventional magnetic systems or spin glasses . 
Most results are for classical systems. The first frus- 
trated system, for which the richness of classical ground 
states was noted, was the triangular lattice d). Subse- 
quently, the pyrochlore lattice, which consists of tetrahe- 
dra that share sites, was identified as a lattice on which 
the frustration effects are especially strong Unusual 
low-energy properties - in particular the absence of order- 
ing at any temperature, was predicted both for discrete 
1^ and continuous |^ classical spin sytems. The ground 
state and low energy properties of the classical pyrochlore 
antiferromagnet - whose quantum version is focus of this 
letter - has been extensivly studied in Q. Both the in- 
terest and difficulty in studying frustrated spin systems 
stem from the large ground-state degeneracy, which pre- 
cludes most perturbative approaches. 

As is the case for most other strong interacting sys- 
tems in more than one-dimension, very little is known 
exactly about the ground states of frustrated quantum 
spin systems. Most of the present knowledge has been 
obtained by numerics or clever approximations. Quan- 
tum fluctuations have been studied in the limts of large- 
S § , where a tendency towards lifting the ground-state 
degeneracy in favor of an ordered state ("quantum or- 
der by disorder") was detected. In the opposite limit - 
S = 1/2, where quantum fluctuations are much stronger 
- the pyrochlore afm has been identified as a candidate 
for a quantum disordered magnet ("quantum spin liq- 
uid") 0, and it has also been discussed in terms of an 
RVB approach However, there are no exact results 
against which to test the reliability of the results in this 
limit. In contrast to this, for conventional ~ bipartite - 
antiferromagnetic spin systems it is well known, for ex- 
ample, that the energy levels are ordered in a natural 
way according to spin, starting from spin zero Geo- 
metrically frustrated systems are not bipartite and thus 
this otherwise quite general theorem does not apply. 

In this letter, we find some first exact results for the 
fully frustrated quantum antiferromagnet on a pyrochlore 
checkerboard, with the help of the reflection symmetry 
of this two dimensional lattice. We shall establish rigor- 



ously that there is always a ground state with total spin 
zero. Furthermore, in the periodic case all ground states 
(if there is more than one) have total spin zero, the spin- 
expectation vanishes for each frustrated unit (this is the 
quantum analog of the "ice ru/e" for the Ising system), 
and there are concrete upper bounds on the susceptibil- 
ity. 

Geometric frustration occurs typically for spin systems 
with interactions that favor anti-alignment and involve 
fully connected units of three or more spins, that can 
obviously not all be mutually anti-aligned. The kagome 
lattice is an example of frustrated spin systems with site- 
sharing triangular units, the pyrochlore lattice and its 
two-dimensional version, the pyrochlore checkerboard are 
examples with site-sharing tetrahedra. We shall focus on 
the latter; it is a two dimensional array of site-sharing 
tetrahedra, whose projection onto a plane is a square lat- 
tice with two extra diagonal bonds on every other square. 
(The regular pyrochlore lattice is a three-dimensional 
array of site-sharing tetrahedra.) The tetrahedra ~ or 
squares with extra diagonal bonds - are the frustrated 
units and will henceforth be called boxes. 




FIG. 1. (a) pyrochlore checkerboard, reflection symmetric 
about dashed line (b) frustrated unit with crossing bonds 

The haniiltonian of a quantum Heisenberg antiferro- 
magnet on a general lattice is (in natural units) 

Hj,^ ^"^Si-Sj, (1) 

where the sum is over bonds that connect sites i 

and j and s = (s^, s^, s^) are spin operators in the spin-s 
representation, where s can be anything. For the checker- 
board lattice the hamiltonian is up to a constant equal 
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to half the sum of the total spin squared of all boxes 
(labelled by x) 

^II(^l+^2+S3+S4)I. (2) 

X 

A 3 X 3 checkerboard with periodic boundary con- 
ditions, i.e., with four independent sites, provides the 
simplest example. It has a hamiltonian that is (up to a 
constant) the total spin squared of one box and the en- 
ergy levels, degeneracies, and eigenstates follow from the 
decompostion of [s]'^'^ into components of total spin; all 
ground states have total spin zero and there are 2s -I- 1 
of them. 

A checkerboard lattice of arbitrary size, with or with- 
out periodic boundary conditions but with an even num- 
ber of independent sites, has the property that it can be 
split into two equal parts that are mirror images of one 
another about a line that cuts bonds, as indicated in fig- 
ure]^, and that contains no sites. We shall now show that 
such a system has at least one spin-zero ground state. It 
is actually not important, for the following argument, 
what the lattice looks like on the left or right; these sub- 
lattices neither need to be checkerboards nor do they have 
to be purely antiferromagnetic (as long as total spin is a 
good quantum number). What is impotant is, that the 
whole system is reflection symmetric about the line that 
separates left and right and that the crossing bonds are of 
checkerboard type. (For a system with periodic bound- 
ary conditions in one direction there will actually be two 
such lines, but we emphasize that PBC is not needed 
here even though it is needed in the usual reflection pos- 
itivity applications; see [|l0| and references therein.) A 
key observation is that these crossing bonds (solid lines 
in figure |l|(b)) form antiferromagnetic bonds Sl • Sr be- 
tween pairs of spins Sl = Si -I- S2 and Sr = 83 + 84 of 
each box on the symmetry line. 

The hamiltonian is H = Hl+Hh + Hc, where Hl and 
Hji act solely in the Hilbert spaces of the left, respectively 
right, subsystem and He contains the crossing bonds. 
For the checkerboard He — J2yii^i + ^2) • (ss + 84))^, 
with the sum over boxes y that are bisected by the sym- 
metry line. Hl and Hfi are completely arbitrary as long 
as they commute with the total spin operator. We will, 
however, assume here that they are real in the ba- 
sis. Any state of the system can be written in terms of a 
matrix c as 

V-^^Cc^jVa "^IV-^)™*, (3) 

where the ^/;^ form a real orthonormal basis of eigen- 
states for the left subsystem and the (?/;^)„t are the cor- 
responding states for the right subsystem, but rotated 
by an angle tt around the 2-direction in spin space. This 
rotation takes t into J,, | into — t; f^nd more generally 
|s,m) into (— )^~™|s, — m). It reverses the signs of the 



operators and S^, while it keeps unchanged. The 
eigenvalue problem Hip = Eip is now a matrix equation 

hLC + c{hnf -EE 4'^c(t«)^ = Ec, (4) 

1=1 y 

where {hi,)ap and {hnjap are real, symmetric matrix ele- 
ments of the corrosponding terms in the hamiltonian and 

(i) 

the ty are the real matrices defined for the spin opera- 
tors Sl and 82 in box y by t^"'^^^ = {ipal^i"'^'' + ^'i'^^'li'f}) 

and t'^^p = i{ipa\s'i'' + S2^^IV'^)- Note the overall minus 
sign of the crossing term in (^): replacing Si -I- 82 by 
83 + 84 and i/)^ by {tp^)rot gives a change of sign for di- 
rections 1 and 3, while the i in the definition of t^^'> gives 
the minus sign for direction 2. 

Consider now the energy expectation in terms of c: 

{tPlH]^) =tTCC^hL + tTC^chR-J2^^''^^ty^<ty^y- (5) 

Since H is left-right symmetric and by assumption real, 
we find that for an eigenstate of H with coefficient ma- 
trix c, there is also an eigenstate with matrix c' and, by 
linearity, with c + c' and i(c — c^). Without loss of gener- 
ality we may, therefore, write eigenstates of H in terms 
of Hermitean c = . We shall also take ip to be normal- 
ized: ('0|'0) = tr c^c = 1. Following let us write the 
trace in the last term of (H) in the diagonal basis of c: 

-trct4'^c(4'V = - Efc,; CfeC/|(4'^)fc;p. This expression 
clearly does not increase if we replace all the by their 
absolute values \ck\, i.e., if we replace the matrix c by 
the positive semidefinite matrix |c| = \/(?. The first two 
terms in (^ and the norm of ?/; remain unchanged under 
this operation. We conclude that if c is a ground state 
than so is |c|. Since c = c"*" — and |c| = + c~ , with 
p.s.d. c+ and c^, we may, in fact, chose a basis of ground 
states with p.s.d. coefficient matrices. 

Next, we will show that the state ■(/'o with the unit 
matrix as coefficient matrix (in the eigenbasis) has 
total spin zero. Since the overlap of tpQ with a state with 
matrix c is simply the trace of c, which is neccessarily 
non-zero for states with a p.s.d. matrix, and because spin 
is a good quantum number of the problem, this will imply 
that there is a least one ground state with total spin zero. 
First consider a spin-1/2 system. In the eigenbasis 
every site has then either spin up or down. The state 
with unit coefficient matrix is a tensor product of singlets 
on corresponding pairs of sites i E L, i' E B. oi the two 
sublattices: 

V'o = (g) (T (T).ot+ i ilU)u' - (g) (Ti - it),,' . (6) 

The analogous state for a system with arbitrary spins, 

s 

^o = (g) E ((-)^""k,m)®k,-m))^^,, (7) 

i^L m— — s 



2 



is also a tensor product of spin-zero states. 

Finally, we would like to show that the projection 
onto the spin zero part of a state with p.s.d. coeffi- 
cient matrix preserves its positivity. This is only of 
academic interest here, but it is non-trivial and may 
very well be important for other physical questions. To 
find the projection onto spin zero we need to decom- 
pose the whole Hilbert space into tensor products of 
the spin components [j]k (8 [j']k' of the two subsystems; 
here fc, k' are additional quantum numbers that distin- 
guish multiple multiplets with the same spin j. Only 
tensor products with j = j' can have a spin zero sub- 
space, which is unique, in fact, and generated by the 



crossing terms of box z in the diagonal basis of c is now 



spin zero state 



m=-j 



\j,m,k)<g>{-y-"'\j,-m,k'). Not- 



ing that (— — m. A;') is precisely the spin-rotated 
state TO, A:'))„t, we convince ourselves that the pro- 
jection onto spin zero amounts to a partial trace over 
TO in a suitably parametrized matrix c. This operation 
preserves positive semidefiniteness, so we actually proved 
that the checkerboard has at least one ground state that 
has both total spin zero and a p.s.d. coefficient matrix c. 

We do not know how many ground states there are. 
To determine the spin of any remaining ground states we 
add an external field to the hamiltonian and study the re- 
sulting magnetization. We will see that the spontaneous 
magnetization of every box on the symmetry line van- 
ishes for all ground states, and thus if we have periodic 
boundary conditions in at least one direction, the total 
magnetization vanishes. Since S'f^^ is a good quantum 
number and S^^ commute with the hamiltonian, this will 
imply that all ground states in such a system have total 
spin zero. Let us thus modify the original hamiltonian 
by replacing the term (s^'^^ -I- s^^ + s^^ + 54^^)2 for a 
single box, z, on the symmetry line by (s^^^* +•§2'^'' -1-33^'' -I- 

si^'' — &)z, i.e., effectively adding a field b to the spins in 
box z and a constant term ^b^ to the hamiltonian. We 

want to distribute the resulting &-terms (s^^'-l-Sj'^^ — 6/2)^, 
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'b/2y, and2(si' 
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to Hl, Hii, He respectively. We cannot use the spin ro- 
tation as before, because the crossing terms in the hamil- 
tonian would no longer be left-right symmetric in the ba- 
sis (|3|). To avoid this problem we will, instead, expand 
eigenstates ip in the same basis on the left and on the 
right: 



(8) 
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In this basis the hamiltonian is left-right symmetric and 
we may assume, as before, that c — . Except for the 
presence of b in box z the energy expectations on the left 
and right are as before. The energy expectation of the 
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This expression clearly does not increase if we replace the 
Ck by their absolute value |cfc| and change the signs of the 
first and last terms. The sign change can be achieved by 
simultaneously performing a spin rotation and changing 
the sign of the field b in the right subsystem. This ac- 
tually completly removes b from the hamiltonian. We 
have thus shown that the ground state energies of the 
systems Hb with and -ffo without the &-terms satisfy 
the inequality Ei, > Eq. Let |6) be a ground state of 
Hi, and |0) a ground state of Hq. It follows from the 
variational principle, that (0|7?{,|0) > {b\Hi,\b) = E^ > 
Eq. Expressed in terms of spin operators this reads 
EQ-2{0\b{sf^+s'i'>+s^^^+sf'>),\0)+b^ >Eq. Recalling 
that we are free to choose both the sign and the magni- 
tude of b wc find that the ground state magnetization of 
box z must be zero: 



(0|(4' 



(3) I J3) , J3) . (3)n 



.10) =0. 



(9) 



This quantum analog of the "ice rule" is true for any box 
on the symmetry line and it holds for all three spin com- 
ponents. In a system with periodic boundary conditions 
and an even number of sites in at least one direction we 
can choose the symmetry line(s) to intersect any given 
box, so in such a system the magnetization is zero both 
for every single box separately and also for the whole 
system: (OjS'tot |0) =0. As mentioned previously this im- 
plies that the total spin is zero for all ground states of 
such a system. 

Let us return to the inequality Ei, > Eq. It implies 
a bound on the local susceptibility of the system: Let 
E{b) = {b\Ho — bSl^i\b) be the ground state energy of 
the periodic pyrochlore checkerboard with a single box 
immersed in an external field b. Recalling Hi, = Hq — 
bSi^l + ^b^, we see that the above inequality becomes 
E{b) + ^b'^ > E{Q) and, assuming differentiability, implies 
an upper bound on the susceptibility at zero field for 
single-box magnetization 



1 d^E{b) 
A db^ 



6=0 
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^4' 



(10) 



where A = 4 is the number of spins in a box. (The sus- 
ceptibility is given in natural units in which we have ab- 
sorbed the g-factor and Bohr magneton in the definition 
of the field b.) 

We would like to get more detailed information about 
the response of the spin system to a global field {b^} in a 
hamiltonian H^i,^ y which is identical to except for the 
terms for the third spin component, which are replaced 



by (4^^ + 4'^ 



,(3) 



,(3) 



bx)^- From what we have 
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seen so far, it is apparent that the corresponding ground 
state energy -E{b^} is extremal for bx = 0. With the 
help of a more sophisticated trace inequality that 
becomes relevant whenever the matrix c in cannot be 
diagonalized, one can actually show that E^i,^-^ has an 



absolute minimum at br = 0: 



£^{b^} > Eq. 



(11) 



Note that we had to put the field on the boxes for this 
result to hold; not every field on the individual spins can 
be written this way. The special choice bx — B/2 cor- 
responds to a global homogenous field B on all spins. 
(The factor 1/2 adjusts for the fact that every spin is 
shared by two boxes.) If E{B) = {B\Ha - BSi!}\B) 
is the ground state energy of the periodic pyrochlore 
checkerboard in the external field B, then (11) implies 



E{B) + ^B^ > E{0), and thus an upper bound on the 
susceptibility per site at zero field (in natural units) 



X A 



1 d^E{B) 
A 



dB^ 



< 



(12) 



where A is the number of independent sites, which equals 
twice the number of boxes. 

All these results continue to hold at finite temperature. 
The analog of ( [Til ) holds also for the partition function 
corresponding to 



Z{b^} < Zq, 



(13) 



as can be shown by a straightforward application of 
lemma 4.1 in section 4 of [|l^ to the pyrochlore checker- 
board. The physically relevant partition function for the 
periodic pyrochlore checkerboard at finite temperature in 
a homogenous external field, Z{B) = tr e~^^-^°~^^tat) ^ 
differs from Zii,\, where bx — B/2, only by a factor cor- 



responding to the constant term in Due to (13) 

the free energy F{B) — —(3^^ In Z{B) satisfies 



F{B) + > ^(0)^ 
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(14) 



This implies (i) that the magnetization at zero field is 
still zero at finite temperature. 



Mt = — r 



1 dF{B) 



A dB 



= 0, 



(15) 
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and, more interestingly, (ii) the same upper bound for 
the susceptibility per site at zero field as we had for the 
ground state: 



XT 



1 d'^F{B) 
A 



dB^ 



< 



1 



(16) 



B=0 



The bounds on the susceptibility hold for arbitrary in- 
trinsic spin-s and agree very well with the results of p^ ] 



for the classical pyrochlore antiferromagnet in the un- 
diluted case. 

It is not essential for our method that only every other 
square of the pyrochlore checkerboard is a frustrated unit, 
only the refiection symmmetry and the antiferromagnetic 
crossing bonds are important. We could, e.g., have di- 
agonal bonds on every square, but then the horizon- 
tal and/or vertical bonds must have twice the coupling 
strength. Our results also apply to various 3-dimensional 
cubic versions of the checkerboard, e.g., with diagonal 
crossing bonds in every other cube. While the method 
does not directly work for the 3D pyrochlore lattice be- 
cause its geometry is too complicated, it has been seen 
in Q that classically this system has similar properties 
to the pyrochlore checkerboard, which is also fully frus- 
trated, and has the added advantage of being more easily 
visualizeable. 

We would like to thank Roderich Moessner for bring- 
ing this problem to our attention, for explaining his work 
on the classical pyrochlore antiferromagnet, and for nu- 
merous helpful discussions and suggestions. 
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